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Abstract 

Multibeam satellite systems have been employed to provide interactive broadband services to geo- 
graphical areas under-served by terrestrial infrastructure. In this context, this paper studies joint multiuser 
linear precoding design in the forward link of fixed multibeam satellite systems. We provide a generic 
optimization framework for linear precoding design to handle any objective functions of data rate with 
general linear and nonlinear power constraints. To achieve this, an iterative algorithm which optimizes 
the precoding vectors and power allocation alternatingly is proposed and most importantly, the proposed 
algorithm is proved to always converge. The proposed optimization algorithm is also applicable to 
nonlinear dirty paper coding. As a special case, a more efficient algorithm is devised to find the 
optimal solution to the problem of maximizing the proportional faimess among served users. In addition, 
the aforementioned problems and algorithms are extended to the case that each terminal has multiple 
co-polarization or dual-polarization antennas. Simulation results demonstrate substantial performance 
improvement of the proposed schemes over conventional multibeam satellite systems, zero-forcing and 
regularized zero-forcing precoding schemes in terms of meeting the traffic demand. The performance of 
the proposed linear precoding scheme is also shown to be very close to the dirty paper coding. 
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1. Introduction 

Multibeam satellite systems have been inspired by the success of the cellular paradigm, which allows 
carefully planned frequency reuse while keeping intercell interference within acceptable limits to achieve 
high spectral efficiency. In addition, the demand for interactive data services on top of broadcasting has 
supported the implementation of multibeam systems, which allow for finer partitioning of the coverage 
area and independent stream transmission within each beam. 

A large number of spotbeams can be employed to cover the same coverage area contrary to recent 
satellite technology where a single (global) beam is employed. Currently, tens or hundreds of beams are 
possible with a typical reuse factor of four. However, due to the antenna design, the beam patterns partially 
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Fig. 1. A multibeam satellite system with a Ground Station serving a 7-beam cluster. Top: Conventional 4-color frequency 
reuse scheme. Bottom: Multibeam joint preceding paradigm using full frequency reuse. 

overlap on the ground creating interbeam interference. The beam patterns and the corresponding allocated 
power have to be carefully designed to ensure that interbeam interference stays within acceptable limits, 
which are determined by the carrier to interference ratio of the beam-edge users. A similar effect has 
been limiting the performance of terrestrial cellular networks for decades, but has been alleviated based 
on multicell joint processing, where user signals in the downlink channel are jointly precoded before 
being transmitted by neighboring BS antennas in order to mitigate inter-cell interference. However, one 
of the practical obstacles in terrestrial implementation is the requirement of a backhaul network which 
enables the cooperation amongst neighboring BSs. 

The principle of multibeam joint processing can be applied to multibeam satellite systems. As illustrated 
in Fig. [T] instead of being served by only one beam, each user's signal is precoded at the gateway (GW) 
and sent by all beams. The main implementation advantage over terrestrial wireless systems is that 
usually the signals for adjacent beams are transmitted from the same GW through the satellite to the 
users in the forward link (FL), as a result, joint precoding can take place at that GW and there is no 
need for expensive backhauling. When multiple GWs serve clusters of beams, distributed joint precoding 
techniques can be employed, but in this paper we focus on one cluster of beams served by a single 
GW. To mitigate interference among multibeams, spatial processing and specifically effective precoding 
techniques can be exploited, which jointly pre-processes data to all beams at the GW. It was reported 
in |[T| that with simple linear precoding, an improvement of the achievable spectral efficiency of about 
25 ~ 50% can be achieved. 

In the following, we provide a review on the available multiuser joint precoding techniques in terrestrial 
wireless communications, as well as an overview of related work in the satellite literature. 

A. Precoding Techniques and Satellite Literature 

In the multiple-input multiple-output (MIMO) and multiple-input single-output (MISO) FL literature, 
a number of linear and non-linear techniques have been proposed to effectively manage interference 
among users. Dirty paper coding (DPC) is an optimal non-linear technique based on known interference 
pre-cancellation which has been shown to achieve the MIMO downlink sum capacity ||2| |[3|. Although 
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DPC is the optimal capacity-achieving transmission strategy, it is not suitable for practical implemen- 
tation and mostly serves as a performance upper bound for the evaluation of less complex techniques. 
Tomlinson-Harasima precoding (THP) is another nonlinear precoding technique which is a more practical 
implementation of DPC based on modulo operations over the constellation symbols |4|. To further reduce 
the implementation complexity, linear precoding schemes have received much attention. Among them, 
zero-forcing (ZF) is one of the simplest techniques that prefilters the transmit signal vector with the 
channel pseudo inverse |[5| so that there exists no inter- user interference. Regularized ZF (R-ZF) [6] 
extends ZF by taking into account the noise variance in order to improve performance in the low SNR 
regime. Opportunistic beamforming (OB) is another linear technique, where each user selects amongst 
predefined random beamformers from a codebook based on both instantaneous channel state information 
(CSI) and long-term channel statistics |[7|. 

In the literature of satellite communications, the majority of early work focuses on the conventional 
scenario where polarization or partial frequency reuse is employed to mitigate interbeam interference. For 
example, authors in |8| optimize the power and beam allocation in order to adapt to channel conditions 
and meet traffic demands. Recently, multibeam joint processing scenarios have been studied in various 
settings. To be specific, FL cases have been investigated in f\], fT^-fl3\, while reverse link cases in 
|[T4)-| 16 1. In all FL studies, fixed satellite services were considered since reliable CSI feedback can only 



be acquired for slow-fading channels due to the long propagation delays. Various characteristics of the 
multibeam satellite channel were taken into account such as beam gain |12|, |13|, rain fading |[T|, 



interference matrix |12| and correlated attenuation areas |13|. In terms of precoding techniques, THP 
was studied in ||9| and p2| , while linear precoding such as ZF and R-ZF were evaluated in |[l|, |[9|. 
Finally, authors in (13] have considered an OB technique based on a codebook of orthonormal precoders 



and low-rate feedback. 



B. Contributions 



In previous work, optimized linear precoding design with concrete objectives and constraints has not 
been treated therefore it is hard to evaluate the potential of each precoding scheme. This paper aims to 
fill this gap and studies the linear precoding design for the FL of a fixed multibeam satellite system with 
perfect CSI using optimization techniques. We take into account the power flexibility, which is essential 
for optimum resource allocation in multibeam satellite communications. It can be implemented by using 



traveling wave tube amplifiers (TWTAs) |22| based on the interaction between an electron beam and the 



radio wave within a tube and the multi-port amplifiers (MPAs) |23|, where the total available power of 
a set of amplifiers can be flexibly distributed amid different beams. We first assume the user terminals 
employ single polarization then extend it to the case of dual polarization, which correspond to multiuser 
MISO and multiuser MIMO systems, respectively. 

The design objective is to optimize any given function of individual users' rates subject to general 
linear and nonlinear power constraints which can model the power flexibility. Specifically, this paper 
makes the following contributions: 
1) We formulate a general precoding design problem by incorporating any functions of individual users' 
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rate and general linear and nonlinear power constraints 

2) We propose a generic iterative algorithm based on alternating optimization to tackle the above 
problem and most importantly prove its convergence. We also devise an efficient algorithm for 
a sub-problem of power minimization. The proposed algorithm is applicable to nonlinear DPC 
precoding. 

3) We extend the proposed algorithm to cope with the case where satellite terminals have multiple 
receive dimensions, and one particular important scenario is when two co-polarization or cross- 
polarization antennas are used at receivers. Polarization correlation and discrimination are considered 
in the channel modeling. 

The main findings via simulations are summarized below 

i) The performance of the optimized linear precoding is very close to that of the optimized nonlinear 
DPC precoding, which provides a performance upper bound; 

ii) To improve the performance, the flexible transmit power enabled by TWTAs and MPAs is even 
more important than the nonlinear precoding. 

The optimization indeed adds to the complexity of signal processing, but this is not a major issue 
since the precoding design is optimized at the GW and thus no additional computation is needed on the 
satellite. 

The remainder of this paper is structured as follows. In section II, we introduce the fixed multibeam 
satellite system model and the problem of linear precoding design with general objectives and general 
linear and nonlinear power constraints. In Section III, we devise a generic iterative algorithm by optimizing 
power allocation and precoding vectors alternatingly and prove its convergence. Section IV modifies the 
proposed algorithm to optimize DPC precoding with fixed encoding order and also as a special case, 
we propose an algorithm to find the optimal solution for the fairness maximization problem with convex 
power constraints. In Section V, the proposed algorithm is extended to deal with the case when satellite 
terminals employ co-polarization or cross-polarization. In Section VI, the proposed multibeam precoding 
is compared to conventional multibeam systems and existing precoding techniques through numerical 
simulations and Section VII concludes the paper. 

C. Notations 

Throughout this paper, the following notations will be adopted. | • | denotes the modulus of a complex 
scalar. Vectors and matrices are represented by bold lowercase and uppercase letters, respectively, and 
II • II is the Frobenius norm. The superscript is used to denote the Hermitian transpose of a vector or 
matrix. A ^ means that matrix A is positive semi-definite. A B and A (gi B denote Hadamai^d 
product and Kronecker product of two matrices, respectively. I denotes an identity matrix. I„ and 1„ 
denote an n x n all-one matrix and an n x 1 vector, respectively. E[-] denotes the expected value of a 
random variable. Finally, x ~ CA/'(m, 0) denotes a vector x of complex Gaussian entries with a mean 
vector of m and a covariance matrix of 0. 

'There are mild conditions about general cost functions and nonlinear power constraints which are satisfied by all reasonable 
performance metrics and will be discussed later. 
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II. System Model 

Consider a FL bent-pipe satellite system for transmitting independent streams to multiple fixed terminals 
through multiple beams. A single GW who has perfect CSI manages a cluster of K adjacent beams on 
ground formed by K antenna feeds (single-feed per beam) on board the satellite and full frequency reuse 
among beams is assumed. This system resembles a multiuser MISO downlink in terrestrial communica- 
tions. By employing a time division multiplexed (TDM) scheme, a single user per beam is served for 
each time slot. The uplink of the feeder link is ideal and no inter-cluster interference is considered. We 
also assume the channel experiences slow fading and perfect CSI is available at the GW. This can be 
realized by feedback/training sent from the terminals via return a return channel, which already exists in 
DVB-S2 (241 . 



A. Satellite Channel Model 

One of the main reasons why satellite communications are challenging and different from terrestrial 
communications is due to the satellite channel characteristics, which need to be properly modeled. The 
satellite channel above 10 GHz operating under line-of-sight (LOS) is subjected to various atmospheric 
fading effects originating in the troposphere, which severely degrade system performance and availability 



1 17|. Among them, rain attenuation is the dominant factor and will be taken into account in our modeling. 
In the following we will describe in detail the satellite channel effects including free space loss, rain 
fading and the beam gain pattern. 

1) Free Space Loss (FSL): Due to the earth curvature and the wide satellite coverage, the free space 
loss in each multibeam will not be identical. In order to model this effect, the FSL coefficient of the kth 



multibeam can be written as |18|: 



= (^) vT5(i? « 

where A is the wavelength and d{k) denotes the distance of the fcth beam center from the center of the 
central beam and do ~ 35786 km. 

2) Rain Fading: To model the rain attenuation effect we use the latest empirical model proposed in 
the ITU- Fj^ Recommendation P.618 1 19|. The distribution of the power gain ^ in dB, = 201ogxo(C)' 



is commonly modeled as a lognormal random variable, i.e.. In (Cde) ~ {fi, a), where fi and a depend 
on the location of the receiver, the frequency of operation, polarization and the elevation angle toward 
the satellite. The corresponding K x 1 rain fading coefficients from all antenna feeds towards a single 
terminal antenna are given in the following vector 

h = ^5e-j'^^" (2) 

where denotes a uniformly distributed phase. The phases from all antenna feeds are hard to differentiate 
and assumed to be identical. This is because we consider a LOS environment and the satellite antenna 



feed spacing is not large enough compared with the communication distance |20|. 



^International Telecommunications Union - Radiocommunications Sector. 
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Since rain attenuation is a slow fading process that exhibits spatial correlation over tens of kms, we 
assume that users undergo the same fading when located within the same beam, but independent fading 
among beams. In other words, we assume that each beam comprises a correlated area p3| | [2T| . 

3) Beam Gain: The link gain matrix defines the average signal to interference-plus-noise ratios (SINR) 
of the each user and it mainly depends on the satellite antenna beam pattern and the user position. Define 
one user's position based on the angle 9 between the beam center and the receiver location with respect 
to the satellite and 9^^^ is its 3-dB angle. Then the beam gain is approximated by |9]: 

h{9,k) = ('4^ + m-^\ (3) 



2u 

where u = 2.07123 sin 0/ sin^sdB and Ji and J3 are the first-kind Bessel function of order 1 and 3. The 
coefficient 6max(^) represents the gain at the kth beam centre when FSL is taken into account as given 
by (Q. 

Collecting one user's beam gain coefficients from all transmit antenna into the K x I vector b, the 
overall channel for that user can be expressed as 

h=yWh0b2. (4) 

B. Signal Model 

Assume the data for user k is with unit average power IE[s|] = 1, VA:, where k = 1,2, ■ ■ ■ , K is 
the index of users. The linear precoding vector = ^/pkWk is used to produce a weighted version 
of Sk before transmission, where pk = ||tfc|p is the transmit power and = is the normaUzed 
precoding vector. The purpose of introducing and -w^ is that in some step of the proposed algorithm, 
they are separately optimized. The on board transmitted signal is the superposition of all users' signal, 
i.e., ^kSk- The received signal at user k can be written as 

K 

Vk = h|.tfcSfc + hl ^ tjSj + Uk (5) 

where is the independent and identically distributed (i.i.d.) zero-mean Gaussian random noise. The 
received SINR for user k is 

E,^k Ihlt,f + NoW Ej^kPM^jl' + ^oW 
where A'o is the noise power density and W is the total bandwidth. The achievable Shannon rate is 

rfc = TFlog(l + rfc),Vfc. (7) 

Denote the achievable rate vector as r = [ri, • • • , r;^]. 

We assume the satellite antenna beams are subject to both linear and nonlinear power constraints as 
described below. 
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C. General Linear Power constraints 



Suppose there are L linear power constraints and the Z— th constraint is expressed as 

K 



(8) 



where (7; > is the power Umit and Q/ ^ is a shaping matrix which includes the following power 
constraints as special cases (we omit the constraint index / for simplicity): 

• Total beam power constraint: Q = I, L = 1; 

• Per-beam (e.g., beam k) power constraint: Qfc = where is a zero matrix except its /c-th 
diagonal element being 1,L = K; 

• Flexible power constraints, e.g., N beam antenna feeds in the set Q = {fci, • • • ^k^} have power 
sharing constraint: Q is a zero matrix except its diagonal elements with indices in Q being 1. 

D. General Nonlinear Power constraints 

Although linear beam power constraints are widely adopted, nonlinear power constraints are also 
relevant. For example, the radio frequency (RF) output power of the high power amplifier (HPA) on 
board the satellite is a nonlinear function of the available direct current (DC) power, which is very costly. 
So if there is a limit on the total DC power, a constraint on nonlinear functions of beam powers is needed. 

Suppose the nonlinear functions linking the beam output power Zk to input power is Zk = gk{xk) 
where gu is a nonlinear function. Here we assume that gu is a continuous and increasing function and 
represents a one-to-one mapping between input and output, thus its inverse function g''^{-) exists and 
Xk = g~^^{zk)- We assume there are J nonlinear power constraints with preset power limits {Pj} on the 
beam antenna output power t|.Dfctfc,VA;, i.e.. 



E. Objective Functions 

Suppose the required traffic demand for user k is F^. We consider a general objective function /(r) 
of rates to optimize. We don't impose any special structure on /(r) but we require that the function is 
continuous, which are satisfied by all commonly used performance metrics. Typical choices include 

1) Throughput maximization: max^^^ r^, which is a widely used performance metric in terrestrial 
communications; 

2) Rate balancing: maxmiiiyt which is another way to meet the traffic demand while maximizing 

^ k 

the worst user's rate to provide proportional fairness. 

3) Rate matching: min^^^ \Fk — r^l", where n is a predefined order. This is especially useful for 
the design of satellite communications |8J to evaluate how well the traffic demand is met. 

Different weights can also be introduced to indicate the priority or delay constraints for each user. We 
impose the additional constraint < F^, i.e., not to over satisfy any users' demand in order to save 

^it may not necessarily be the input power to tlie amplifier and could be any power of interest 



K 




(9) 



k=l 
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costly on board power and this has been included automatically in the rate matching problem. With this 
additional constraint, it is easily seen that throughput maximization is a special case of rate matching 
when n = 1 as 

K K K 

Y.\^k-rk\=Y,Fk-Y,rk. (10) 

k=l k=l k=l 

Note that different objectives in general result in different optimized rate vectors and the choice of the 
appropriate objective is up to the satellite operator. 

F. Problem Formulation 

The problem of interest is to minimize a general objective function /(r) of rates subjective to general 
linear and nonlinear power constraints by designing precoding vectors {tfc,VA;}. Mathematically, it is 
expressed as 

min fir) (11) 
K 



s.t. Y^t\ciitk<qi,yi, 



k=l 
K 

T.9kjiAWk)<Pj,yj, 

k=l 

Tk < Fk,yk. 

The main difficulty lies in the fact that the objective function /(r) is non-convex in general, e.g., the 
throughput maximization problem is non-convex and has been studied in terrestrial communications p5|. 
The nonlinear (possibly non-convex) power constraints also introduces additional difficulty. 

III. Generic Precoding Design For Rate Enhancements 

To tackle the difficulty of non-convex objective function and nonlinear power constraints, we adopt 
the strategy of alternating optimization [25], that is, we first optimize the precoding vectors {w^} then 
update the power allocation {pk} with fixed precoding vectors. 

A. Optimization of Precoding Vectors {w^} 

Optimizing the precoding vectors under the non-convex objective function and non-linear power 
constraints is too complex. Instead we propose to solve the following per-beam power minimization 
problem subject to minimum rate constraints, general linear power constraints (first kind) and per beam 
power constraints which are also linear (second kind): 

min 7 (12) 
S.t. rk>Rk,yk, (13) 

K 

^t^.Q/tfc < qi,yi, 
k=l 
K 

Y,tinjtk<lPj,j = !,■■■ ,K, (14) 

A;=l 
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where Pj is a known parameter related to the j-th beam's power and 7 is an auxiliary variable. The 
above problem aims to reduce each beam's transmit power while satisfying the rate constraints {Rk}, 
which is a predefined minimum rate requirement. The reason of such a formulation will become clear 
as our discussion goes on. 

Note that the above optimization problem is convex since the objective function is convex, the two 



kinds of linear power constraints are convex and the individual rate constraints ( 13l, although originally 



non-convex, can be made convex using the transformation in |I261| and rewritten as 



(2^ 




t|hfc|2 + WNo 



<htt 



(15) 



and as a result problem (12 1 is recognized as a second-order cone programming (SOCP) problem which 
is convex and readily solved using standard techniques |27|. 

An interesting fact about (12i is that at the optimum, the rate constraints Rk = rj.^'ik, which aligns 



automatically with the last constraint in the original problem (111 



Once we get the optimal t^, the normalized precoding vector can be obtained as v^k = p^- The power 
is also optimized in ([12]) via t^, although it does not take into account the general objective function 
/(r), it is still very useful in developing the algorithm and for this purpose the beam power is stored as 
Pk = ||wfc|p,V/c for later use. 



Notice that the problem (12) is also useful to detect whether the traffic demand can be satisfied or 



not, i.e., replace R^ with Fj. and solve ( 12 1. If it is feasible, then the optimization is done as the traffic 
demand is perfectly met. In the remainder of this paper, we always assume that the traffic demand is 
high enough such that it is not possible to satisfy all users' rate constraints simultaneously. 



B. Power Optimization 

With fixed precoding vectors {w^}, the general power optimization problem can be modified from 



(111 and formulated as 



mm 

{p4 



s.t. 



/(r) 



K 

k=l 
K 

E5fcj(4Ditfc)<p„vj, 
k=i 



(16) 



rk = Wlog2 1 + 



< Fk,yk. 



E,^kPj\K^j\' + NoW^ 

It is seen that both the linear power constraints and the maximum rate constraints are in linear form 
about power {pk} which are easy to handle. However, due to the possibly non-convex objective function 



/(r) and nonlinear non-convex power constraints, ( 16) is non-convex in general and its globally optimal 
solution is difficult to find. 
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Here we propose to find its local optimal solution by using any gradient-based numerical algorithm, such 
as the steepest descent algorithm, sequential quadratic programming method and trust region optimization 



28|. Although the globally optimal power solution is not guaranteed to (16 1, it will be seen later that 



by properly choosing the initial point, a convergent iteration algorithm can be constructed. 

C. The Proposed Generic Iterative Algorithm and Proof of Convergence 

With the above strategies to optimize precoding vectors and power allocation altematingly, we propose 
a generic iterative algorithm below to jointly optimize the precoding vectors and power allocation. 



Proposed Generic MISO Algorithm 

Step 1: Initialize {w^} and {pk} such that both the linear power constraints and maximum rate 



constraints are satisfied in ( 11 1. Towards this, a convenient choice is that {p^] are chosen 
to be very small values. 
Step 2: For each k, evaluate the achievable rate for user k using ([6]|7]| and store it as R^. 



Step 3: With the above {Rt] as inputs, solve the optimization problem (12i for t^, then obtain 



the precoding vector w^. = the power solution p^ = ||tjt|p, and i\ = t|,Djtfc,V/c. 



Step 4: With fixed {v^k}, solve the optimization problem (11 1 for power allocation {p^} using a 

gradient-based algorithm with {pk} being the initial power solution. 
Step 5: Go back to Step 2 until it converges. 



Although there are many possibilities to construct iterative algorithms based on alternating optimization, 
the convergence is not guaranteed. Next we will prove that the proposed generic MISO algorithm indeed 
converges. 

Theorem 1: The proposed Algorithm 1 always converges. 
Proof: 

We denote the value of objective function /(r) at Step 2,3,4 of the nth iteration as and 
respectively. 

Due to the small initial power {pk}, Rk < Fk,yk always holds true. From Step 2 to Step 3, it is easy 



to see that the problem (12) at Step 3 is always feasible and Rk is satisfied with equality. Therefore we 
have 

(17) 

The improvement at Step 3 is that less or at most equal power is used for each beam as that at Step 2, 
which in turn satisfies all nonlinear beam power constraints (due to their monotonicity) and creates room 



for further power optimization. This is realized by introducing the per beam power constraints ( 14 1. 

At Step 4, the power is further optimized locally and the power solution {pk} at Step 3 is taken as an 
initial point such that the optimized power {p^,} should be no worse than {pk} which is obtained based 
on another criterion of minimizing beam power rather than /(r), in other words, 

) < ). (18) 



Combining ( 17 1 and ( 18 1, we have 
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which impUes that the value of the objective function is monotonically decreasing from one iteration to 
the next and as it is lower bounded by zero, the algorithm always converges. This completes the proof. 



We have the following remarks on solving problem (12) at Step 3 of the proposed algorithm, which 
is a key step in the algorithm. 

Remark 1: We adopt the maximum beam power as the objective function in ( [T2] ) and therefore the 
outcome of the algorithm is that we minimize the general objective function /(r) with minimum beam 
power. This choice of objective is general enough while for a specific problem, it is not unique, e.g., if 
there are only linear power constraints, the objective function at Step 3 can be replaced by the maximum 
user power or total transmit power which indicates their respective importance. With these changes of 
objective function at Step 3, the overall convergence is still guaranteed which can be seen from the proof 
of Theorem 1. 



Remark 2: Although (12 1 can be solved using standard numerical algorithm, there are less complex 
algorithms by making use of its analytical structure. We propose a more efficient algorithm in the 
Appendix to solve ( [T2] ) based on its dual problem and subgradient method inspired by [29j . The proposed 
algorithm in the Appendix can be modified to deal with different objective functions. 

IV. An Extension To DPC and A Special Case For Fairness Maximization 

A. Application to DPC 

DPC is known to be the sum capacity-achieving technique in multiuser MIMO downlink. Here we 
will show that given a decoding order, the proposed MISO algorithm can be applied to the precoding 
design in DPC straightforwardly. 

Let us assume that ttq = {1,2,..., K} is a trivial user encoding order. Then the received SINR at user 
k is 

Tfc = ^^^^ (20) 

Z,>kKtj\' + NoW 

which possesses similar structure as (|6]l and thus the proposed algorithm can be easily modified to find 
the precoding vectors for DPC. 

Since determining the optimal encoding order is computationally expensive, we investigate the perfor- 
mance of some intuitive heuristics. In general, users encoded first receive a large amount of interference 
which gradually decreases as the serial encoding proceeds. In this respect, users with good channel 
conditions should be encoded first because they can afford being interfered. Similarly, users with high 
rate requirements should be decoded last because they have to avoid interference to achieve the requested 
rates. As a result, the users are ordered based on the increasing order of the metric: | ^^^ (i+'||hfc|p) }• 

B. Special Case: Optimal and Efficient Solution to Fairness Maximization with Convex Power Constraints 

As mentioned before, in general the objective function /(r) is non-convex which makes the whole 
problem difficult even with linear power constraints only. There has been a special case well studied in 
terrestrial communications |30] where the efficient and globally optimal solution exists: to maximize the 
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worst user's weighted rate in order to achieve fairness subject to the total power constraint. Here we 
extend it to the nonlinear but convex power constraints. Mathematically, the problem is formulated as 

^'^ (21) 



max 
{t4 



mm 

k Fk 

K 



S.t. ^tlciitk<qi,yi, 



k=l 
K 

fc=l 



where we assume g^^i') is a convex function. Next we propose an algorithm within the framework of 
optimization to achieve the optimum. We first introduce an auxiliary variable 7 and reformulate it into 



max 

{tk},i<i 



s.t. 



(22) 



tlhk > 



K 



\ 



(2 



1) J]|t]hfc|2 + TyiVo , vfc, 



k=l 
K 



E5fc](4D,tfc)<p„vj. 



fc=i 



Now ([22]) is a convex problem except for the variable 7. To solve it, we use the bi-section search 
approach [21] to find the optimum solution. To be specific, at each search iteration for a given 7 < 1 



(this guarantees < Fk^Mk), we check whether all constraints in (22 1 are feasible or not. If it is feasible, 
then increase the value of 7 and decrease it otherwise. 



V. Multiple Antennas at Receive Terminals 

In this section, we assume the user terminals have multiple antennas and the system turns out to 
be a multiuser MIMO downlink system. Specifically we study the effects of co-polarization and dual 
polarization. Dual polarization has been proposed as a viable solution for employing single-user MIMO 
techniques in mobile satellite systems | [3T| , but its performance in fixed multibeam satellite systems is 
not well understood. In this section, we investigate the effect of dual-polarization antennas with linear 
precoding in the FL of fixed multibeam satellite systems. 

Assume that each terminal is equipped with multiple antennas. More specifically, the following four 
cases are considered: 1) two co-polarization antennas with receive beamforming (RBF), 2) two cross- 
polarization antennas with RBF, 3) two co-polarization antennas with antenna selection, 4) two cross- 
polarization antennas with polarization selection. In general,RBF is optimal but antenna selection is 
less expensive in computation and implementation. In terms of polarization, cross-polarization antennas 
achieve lower receive SNR because of the power imbalance but they have better de-correlation properties 
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due to propagation characteristics. As a result it is not straightforward to judge which scheme is better 
and the main aim is to choose the best technique based on realistic channel parameters. 

A. Channel Modeling 

1) Two co-polarization antennas with receive beamforming: Assume each of K terminals has two 
receive antennas. The channel for each terminal can be denoted as with dimensions 2 x K. A main 
assumption of the considered channel model is full transmit-side correlation due to the lack of scatterers 
close to the satellite. For the co-polarization case, the two receive antennas are also assumed to be fully 
correlated and thus the channel for the kth user can be written as 



Hit 



I2 (K) hk 



(23) 



where is defined in Q. The beam gain coefficients for each user vary across the beams due to the 
beam pattern, while rain fading is identical due to identical illumination angle. 

2) Two cross-polarization antennas with receive beamforming: Assume that both satellite feeds and 
terminal antennas are dual-polarized. For the sake of simplicity, a Kronecker correlation model between 
polarizations is assumed with both transmit and receive correlation. In this direction, the channel for the 
kth user can be written as the 2 x 2N matrix: 



Hi 



\/Vax(fc) ih^H) 



R^Hz^Ri 



1 



N 



(24) 



where P is a 2 x 2 matrix modelling the power imbalance between polarizations, H^ includes the rain 
fading coefficients for co- and cross-polarization subchannels and R is the correlation matrix between 
polarization^ More specifically, the matrix Hjt can be written as: 



Ht 



(25) 



■^00 ^01 

A A. 

where = {0, 1} are assumed to be independent rain fading coefficients as modelled in Section 

II.A. (/) is a uniformly distributed phase. Furthermore, if a is the power imbalance factor and p is the 
correlation factor, P and R can be modelled as: 

(26) 



R 



1 


a 


a 


1 


1 


P 


P 


1 



(27) 



It should be noted that the a factor is related to Cross Polarization Discrimination (XPD) as follows: 



''it should be noted that due to the statistics of the rain fading coefficients (non-zero mean), the channel power of will 
be higher than the one of Hfc. In order to ensure a fair comparison amongst all scenarios, the channel matrices are normalized 



so that trace(HfcHj.) = trace(HfcH 
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B. Receive Strategy I: Antenna Selection 

1) Two co-polarization antennas with antenna selection: Since in the co-polarization case the two 
antennas are considered fully correlated, this scenario degrades to the single-antenna terminal already 
studied in the previous sections. 

2) Two cross-polarization antennas with polarization selection: Building on scenario 2 we discern two 
cases. In the first case, each satellite antenna has single polarization. The corresponding 2 x N channel 
matrix can be easily derived from ( [24] ) by picking alternating polarizations. In the second case, we also 
select the best polarization at the terminal using the criterion 

argmax ||Hfc(n, :)||, V/c (28) 

n 

which results in a 1 x A'^ channel matrix. 

Both of the above cases result in equivalent multiuser MISO channels, therefore the proposed generic 
MISO algorithm could be applied to find the optimized precoding vectors. 



C. Receive Strategy II: RBF and The Proposed Generic MIMO Algorithm 

When there are multiple receive antennas at the terminal, assuming single-user decoding, the optimal 
strategy is to use RBF vector {u^} that matches the channels and the precoders. The received SINR for 
user k can be expressed as 



|u|.Hfctfcp 



(29) 



and for fixed precoding vectors t, the optimal RBF vector to maximize the received SINR is given by 



K 



-1 



(30) 



With the above result, we propose the following generic MIMO algorithm based on the proposed 
ge neric MISO algorithm. 



Proposed Generic MIMO Algorithm 

Step 1: Initialize {w^fc,Ufc,pfc} such that both the linear power constraints and maximum rate 



constraints are satisfied in (11 1. To guarantee that, normally {pk} need to be very small 
values. 



Step 2: 
Step 3: 



Step 4: 
Step 5: 
Step 6: 
Step 7: 



For each k, evaluate the achievable rate for user k using (29 1 and store it as Rk- 

With the above {Rk] as inputs and given fixed {u^}, solve the optimization problem 



similar to (12) with substitution h 



r, the power solution 



u^HfcjVA; for {t^}, then obtain the precoding 

ti 



and Pk = t[Djtfc,V/c. 



vector = 
Update the RBF vector {ufc} using (30 1. 

Evaluate the achievable rate for user k using (29 1 and store it as Rj^,\/k. 
If Rk < Fk, then update u^. = u^^Mk. 



With fixed {wfc,Ufc}, solve the optimization problem similar to (111 with substitution 



= u|,Hfc,V/c for power allocation {p^} using a gradient-based algorithm with {p^} 



being the initial power solution. 
Step 8: Go back to Step 2 until it converges. 
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The difference between the proposed MIMO algorithm and the proposed MISO algorithm is due to 
Step 4-6, which update the RBF vectors and ensures the rate requirements are not over satisfied. The 
proposed generic MIMO Algorithm can be proved to always converge using the same argument in the 
proof of Theorem 1 and thus the proof is omitted here. 



VI. Simulation Results 
Computer simulations are conducted to evaluate the performance of the proposed algorithms. For 



multibeam joint precoding, we consider a fixed satellite system as described in Section II-A with detailed 
parameters listed in Table |l] The rain fading corresponds to temperate central European climate. We 
assume there are K=7 on board antenna feeds serving 7 beams on the ground. Within each beam, there 
are 4 fixed user terminals and they are served in a TDM manner. Achievable rates for all users will 
be shown as the performance metric for different schemes and different objectives. Unless otherwise 
specified, we choose the I2 norm minimization, /(r) = J2k=i — ''fcP. as the objective function. 
Traffic demand is assumed to be asymmetric and uniformly distributed with the mean vector listed in 
Table I. For linear power constraints, as listed in Table |l| we assume each satellite RF saturation power 
is 80 dBW, which is the individual beam power constraint. 



TABLE I 

Satellite Scenario Parameters 



Parameter 


Value 


Orbit 


GEO 


Frequency band 


20 GHz 


Number of beams 


K = 7 


Beam diameter 


D = 250 km 


3dB angle 


eadB = 0.4° 


Rain fading mean 


At = -2.6 


Rain fading variance 


cr = L63 


Polarization 


Single/Dual 


Max antenna Tx gain 


52dBi 


TWTA RF power @ saturation 


SOW 


User terminal maximum antenna gain 


41.7dBi 


FL free space loss 


210dB 


User link bandwidth 


W=500MHz 


Clear sky receiver temperature 


207°K 


Mean of uniformly distributed traffic demand 


[4 0.8 0.8 0.8 2 2 2] Gbps 


Frequency reuse factor for the conventional scheme 


4 



For the MISO case where the user terminal has only one receive antenna and single polarization, the 
proposed generic MISO optimization algorithm will be compared with the following schemes: 
1) Conventional schemes with single-beam processing, the same bandwidth W and frequency reuse 
factor 4. The achievable rate for user k is 
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where pk is the transmit power for beam k and co-channel beam interference is ignored; 

2) ZF preceding: collect all users's channels into a K x K matrix H = [h|; • • • ; hj^], then preceding 
vector Wk is taken from the normalized A;-th column of (H1^H)~^H1^ 

3) R-ZF preceding: Wk is taken from the normalized k-th column of (Hl^H + aI)~^Hl^ where a = 
1 33 1 and Pq = 80 W is the maximum beam power constraint; 

4) DPC with the encoding order determined by | (i+'||hfc||^) } nonlinear precoders are 
optimized using the proposed algorithms to provide a performance upper bound. 

Performance results are depicted for all 7 beams whose indices are shown on x-axis. Fig. |2] shows 
individual users' rates for different MISO schemes to match the traffic demand. The throughput and the 
average value of the I2 norm cost function (n = 2 in the rate matching function in Section lI.E) are shown 
in the legend. First it is verified that all multibeam schemes outperform the conventional single-beam 
processing. The simple ZF preceding can achieve 67% more throughput than the conventional scheme. 
R-ZF performs slightly better than ZF by considering the noise effects. Compared with ZF precoders, it is 
observed that the proposed generic optimization scheme achieves slightly higher throughput but matches 
the traffic much better (40% lower I2 norm). It is also seen that DPC only has very marginal performance 
gain (around 4% in terms of throughput and 6% for I2 norm) over the proposed linear preceding. We 
also study the impact of flexible total power constraint of 7 x 80 = 560 W for the proposed scheme and 
DPC. It is seen that higher rates are achieved compared with the individual power constraints, which is 
expected. The proposed scheme with the total power constraint even outperforms DPC with individual 
power constraints, which implies that flexible power constraints are more important than the complicated 
non-linear signal processing. It is again observed that the proposed linear preceding performs nearly as 
well as DPC. This indicates that linear processing is adequate for multibeam satellite. 

Fig. [3] shows the average actual power used by each beam for different transmission schemes. The 
total power is also shown in the legend. First it is noted that all schemes have a tendency to allocate 
power adaptively to traffic demand. The conventional scheme always uses the highest amount of power 
due to the lack of multibeam cooperation. ZF solution uses only about 60% total power of that used by 
the conventional scheme which is a huge saving of satellite power. The proposed scheme consumes less 
power but more adaptively than ZF to meet the traffic demand within the per beam power limit. To get a 
clearer picture of achieved rate and used power. Fig. |4] shows the rate efficiency normalized by the actual 
power used by each beam. It is confirmed that the conventional scheme has a very low rate efficiency 
due the low rate achieved and high power consumed. ZF solution has a high efficiency for those beams 
with low traffic demand due to the low power used. 

Fig. [5] illustrates the impacts of choosing rate balancing, I2 norm minimization, throughput maximiza- 
tion as objective functions. It is seen that although the objective of rate balancing guarantees fairness 
among users, it degrades too much the throughput performance. The proposed schemes with sum rate 
maximization and I2 minimization as objectives greatly outperform rate balancing and achieve 10% and 
8% higher throughput, respectively. Compared with DPC preceding which is optimal for throughput 
maximization, the performance degradation of the proposed linear preceding is almost negligible. 

To evaluate the performance of the MIMO case when user terminals have either co-polarization or 
dual-polarization antennas, the proposed generic MISO and MIMO optimization algorithms are applied 
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to the following cases: 



Case 


Feeds 


Terminals 


Technique 


1 


S ingle-polarization 


2 co-polarization antennas, p = 


1 


RBF 


2 


Dual-polarization 


2 dual-polarization antennas, p = 


0.9 


RBF 


3 


Alternating single-polarization^ 




2 dual-polarization antennas, p = 


0.9 


RBF 


4 


Single -polarization selection 


2 dual-polarization antennas, p = 


0.9 


Best polarization selection 



The results are shown in Fig. [6] The throughput and the averaged value of the I2 norm cost function are 
shown in the legend. It is observed that with co-polarization receive antennas, the I2 norm is significantly 
reduced compared with the MISO case, which indicates better traffic matching. The scheme with dual- 
polarization at satellite antenna feeds and terminals employing RBF achieves the best performance at the 
cost of increased hardware complexity. The performance of the scheme with alternating single-polarization 
satellite and dual-polarization terminals employing RBF is comparable to that with co-polarization receive 
antennas, while polarization selection further degrades the performance. The performance of all dual 
polarization techniques is greatly degraded by the correlation effects. The co-polarization receive antennas 
at the terminal seem to be a more promising solution considering both performance and hardware 
complexity. 

VII. Conclusions 

This paper has studied the optimization of linear precoding for multibeam satellites with general linear 
and nonlinear power constraints. A generic iterative algorithm has been proposed to handle any objective 
function of individual rates and the convergence is proved. The proposed algorithm has been extended 
to design DPC precoding with fixed encoding order and to the case when user terminals are equipped 
with two co-polarization or dual-polarization receive antennas. Simulation results have shown substantial 
performance gain compared with conventional single beam based processing and existing precoders. The 
optimized linear precoding is also shown to be as effective as DPC precoding, which confirms its practical 
value. The impact of co/dual-polarization has also been demonstrated. 

Future work includes the study of of promising non-linear precoding design like THP in a DVB-S2 
system. In addition, in this paper, the feeder link and CSI at the GW are assumed perfect, which is not 
very practical and will be investigated in the future. 

^Each feed is single polarization but we alternate between co- and cross-polarization as we enumerate the feeds. 
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Appendix 



Efficient Algorithm to Solve ( 12 i 



The dual problem of ( 12 1 is 



K 



max 
a>o,A>o,yLt>o 



k=l 1=1 

K L K f I \ 

s.t. ^^^i + ^'Q' + E "^k^^l ^ ^ + — 

i=i /=i k=i ^ 2»' -1/ 



(32) 



K K 



The strategy is to first fix A and fj, and solve (32i for a and find corresponding primary solution {t^} 
then use it to update A and /x. 
Given fixed A, ([32]) reduces to 

K 

max WNq V Uk (33) 
OL>0 ^-^ 



k=l 



K L K / 1 \ 

S.t. Y Z^jDi + Y ^'Q' + E ""^^H ^ 1 + ^khl 

j=i 1=1 ^ 2w -1/ 



V/c. 



whose optimal solution has the interpretation of virtual uplink power and is found by performing the 
following iterative update p2] |: 

1 



oik 



(l + ^) hi (l + Ell + Ek=i «fch,hl) hk 



(34) 



Then the optimal precoding must be in the direction of 

(K L K 

^^jDj + ^ AiQi + ^afchfch^ j hk 



(35) 



1=1 k=l 

and the downlink power can be obtained by using the fact that the rate constraints are satisfied with 



equalities in ( 12 1. 

Summarizing the above results, we propose an efficient algorithm using similar subgradient method in 
1 29 1 as follows. For more details, please refer to |29|. 



19 



Step 


1: 


Step 


2: 


Step 


3: 


Step 


4: 



Proposed Efficient Algorithm to Solve ([T2)l 

Initialize A > 0, /x > 0. 



Given A, /x, find the optimal solution to (33) by using fixed point update d34b iteratively. 



Update tfc using ([35j), V/c. 
Find downlink power using 



where G 



kk 



1 



itihfc|2,G 



kj 



Step 5: 
Step 6: 

Step 7: 

Step 8: 

Step 9: 



2Tr -1 

Update downlink precoding vectors 



|t]hfc|2,Vi/ A;. 



(36) 



Update A; = max(0, A/ + 77/AA;), V/ where rji is the step size. AA/ is the subgradient and 
one possible choice is AA; = Ylik=i ^iQi^fc ~ V/. 

Update = max(0,^j + pjAfij),\/l where pj is the step size. Afij is the subgradient 
and one possible choice is Apj = J2k=i tfeD/t^, j = 1, • • • ,K. 



Find vaUd /x using the projection method into the feasible set {/x : XljLi f^j^j — 

Ef.im- 

Go back to Step 2 until it converges. 
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Fig. 2. Comparison of rates for different scliemes witli tlie objective of I2 norm minimization. In tiie legend, tlie name of eacli 
sciieme is followed by the throughput and the average value of the h norm. 
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Fig. 3. Comparison of beam powers for different schemes. In the legend, the name of each scheme is followed by the total 
beam power. 
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Fig. 4. Comparison of rate normalized by beam power for different schemes. 
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Fig. 5. Comparison of rate with different objectives. In the legend, the name of each scheme is followed by the throughput. 
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Fig. 6. The effects of co-polarization and dual-polarization antennas with correlation coefficient of 0.9 for different schemes. 
In the legend, the name of each scheme is followed by the throughput and the average value of the I2 norm. 



